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An analysis of the P2S6 cluster electronic structure and its comparison with the crystal valence 
band in the paraelectric and ferroelectric phases has been done by first-principles calculations for 
Sn2?2S6 ferroelectrics. The origin of ferroelectricity has been outlined. It was established that 
the spontaneous polarization follows from the stereochemical activity of the electron lone pair of 
tin cations what is determined by hybridization with P2S6 molecular orbitals. The chemical bonds 
covalence increase and rearrangement are related to the valence band changes at transition from the 
paraelectric phase to the ferroelectric one. 
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INTRODUCTION 

For such perovskite ferroelectrics as BaTiOs, the main 
origin of spontaneous polarization is commonly related 
to the hybridization interaction between the transition- 
metal and oxygen ions 1 . Another mechanism involves 
cations with "lone pair" electrons which have a formal 
ns 2 valence electron configuration 2 -. In the same manner 
as for the d° transition-metal ions, these p° ions (as ex- 
ample Pb 2+ for PbTi0 3 , or Bi 3+ for BiMn0 3 ) contain 
some p-charge density which contribute to the displacive 
distortions. If the lowering of energy associated with the 
hybridization interaction is larger than the interionic re- 
pulsion opposing the ion shift, then a ferroelectric distor- 
tion appears. This "stereochemical activity of the lone 
pair" is the driving force for off-center distortion in fer- 
roelectrics. Both named origins of ferroelectricity (the 
first - d -" ness" , and the second - " lone pair" activity) 
are familiar to the second-order Jahn- Teller (SOJT) ef- 
fect*^. This effect is determined by a balance of positive 
and negative contributions to the total energy. The first 
one describes short range repulsive forces and is related 
to the rigid ions (with frozen electronic configuration) 
shifts from original high symmetry positions. Such term 
is small for the cases of "closed-shell" d° or p° cations. 
The second, negative, contribution describes the relax- 
ation of electronic configuration in response to the ions 
displacements through covalent bonds formation. This 
term favors the ferroelectric distortion. For full picture, 
the geometrical (or hybrid improper) mechanism, which 
is related to the rotational modes that trigger instability 
of polar mode 5,6 , could also be considered at study of 
ferroelectricity nature. 

The cubic crystal lattices of ABO3 compounds are 
built by covalent bonds A — O and B — O with con- 
siderable contribution of ionicity for former ones. Nat- 
urally more complex bounding evolution could be sup- 
posed at the ferroelectric phase transition in ion-covalent 
crystal Sn2P2S6 with monoclinic lattice. For this com- 
pound the Sn 2+ cations and the (P2S6) 4- anion clusters 
are joined by mostly ionic Sn — S bonds at covalent 
P — S and P — P bonding. Sn2P2S6 uniaxial ferro- 
electric undergoes the second order phase transition at 



To ~ 337 K (P2i/c — » Pc, two formula units in the 
elementary cell for both phases (Fig. [I|) in a crossover 
displacive-order /disorder region 7 . Ferroelectric instabil- 





FIG. 1. The crystal structure of Sn2?2S6 ferroelectric phased. 
The tin atoms positions in the paraelectric phased are shown 
by red. The shape of the Brillouin zone with denoted sym- 
metrical points is shown for primitive monoclinic lattice. 

ity in Sn2P2S6 crystal is a result of non-linear A Q B\ cou- 
pling of the soft polar B u and fully symmetrical A g optic 
modes, leading to three-well potential^. Here the oppo- 
site picture appears in compare with perovskites where 
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the only one lattice mode could determine dynamical in- 
stability related to the ferroelectric phase transition 10 . 
In general, all the 13 B u and 15 A g optic modes were 
accounted in the frozen phonons approximation for con- 
struction of effective Hamiltonian for Sn2P2S6 crystal 
what was applied for the MC simulation of the ferro- 
electric phase transition and their behavior under hydro- 
static pressure 11 . A system with a three- well potential 
was early considered by Line a 12 i 13 with application for 
the LiNbC>3 an d LiTaC>3 crystals. Such system can be 
generally described by two order parameters (related to 
dipole and quadruple moments), and as a result, a variety 
of stable, metastable and unstable states can be realized 
on a phase diagram o 14 i 15 . 

The strong anharmonicity of Sn2P2S6 crystal lattice 
is obviously joined with effective electron-phonon inter- 
action, that appears as a stereochemical activity of the 
tin cations electron lone pair 5s 2 , and in fact, it is a re- 
flection of the SOJT effect. Possible leading role of the 
cations' stereoactivity for Sn2?2S6 crystals was early 
noted at X-ray structure investigations of their paraelec- 
tric phase in comparison with structure data for the fer- 
roelectric phased. Structural evidences of the tin cation 
stereoactivity was also analyzed in details by structure re- 
finement of the paraelectric and the ferroelectric phases 
for Sn2P2Se6 selenide analog 16 . By Mossbauer effect 
investigations for 119 Sn nucleus^ and by NMR spec- 
troscopy for isotopes 31 P and 119 S n 18 i 19 , the important 
changes of chemical bonding at the ferroelectric phase 
transition in Sn2?2S6 were found. The X-ray photoelec- 
tron spectroscopy confirms growth of the chemical bonds 
covalence in the ferroelectric phase 20 . 

Sn2?2S6 crystals are ferroelectric-semiconductors with 
promising photorefractive 2 ^, photovoltaic* 2 ^, electroop- 
tio^ 3 - and piezoelectric* 2 ^ characteristics. Their ferro- 
electric properties are effectively influenced by state of 
electronic subsystem 2 ^. Influence of the sulfur and 
tin vacancies on semiconductive and optic properties 
of Sn2P2S6 crystals was studied recently 2 ^. These 
data motivate the electronic structure investigation for 
Sn2P2S6 crystals in the paraelectric and the ferroelectric 
phases. 

The first-principles calculations in LDA approach of 
Density Functional Theory (DFT) for Sn2?2S6 ferroelec- 
tric phase were carried out by several groups 2 ^—. By 
Grigas et a/j 20 i 30 , the electronic structure of both para- 
electric and ferroelectric phases of Sn2P2S6 was calcu- 
lated in the cluster approach. For the Sn2P2Se6 selenide 
analog, the electronic structure have been investigated 31 
by first principles calculations only for the paraelectric 
phase. The electronic structure and phonon spectra pres- 
sure dependence for acentric layered rhombohedral crys- 
tal SnP2S6 were investigated theoretically in LDA ap- 
proach 3 - 2 -. For this compound, the tin cations are almost 
fully ionized (Sn 4+ charge state) what exclude possibility 
of stereochemical activity of their 5s 2 electron lone pair. 
The electronic structure of high charged (P2S6) 4- and 
(P2Se6) 4_ anion clusters was discussed in paper a 29 i 31 i 32 



at analysis of Sn2P2S6, Sn2P2Se6 and SnP2S6 electron 
energy spectra. Analysis of P2S6 cluster chemical bound- 
ing was also done in Hartree-Fock approach 33 . These an- 
ions arrangement have been investigated experimentally 
and theoretically in different approximation for the lay- 
ered crystals like M 2 P 2 S 6 (M - Fe, Ni, Mn, . . . 
For the CuInP2Se6 layered compound with two differ- 
ently charged cations, the SOJT effect was established 
as an origin of cooper ferrielectric ordering 39 . 

In this paper, the first-principles calculations in 
LDA approach of DFT for electronic structure of 
Sn2P2S6 crystal in the paraelectric and the ferroelec- 
tric phases were used for analysis of chemical bonds 
transformation at the spontaneous polarization appear- 
ance and for establishing of the ferroelectric state ori- 
gin in the phosphorus containing chalcogenides. As a 
background of the investigations, the free P2S6 structure 
group (atomic cluster) electron spectra and peculiarities 
of their molecular orbitals were considered. The electron 
structure and chemical bonding nature in the paraelectric 
phase will be also analyzed. The stereochemical activity 
of the electron lone pair of tin cations is examined in 
detail. The growth of covalence and ions recharging are 
related to the structure spontaneous polarization. Fi- 
nally, the influence of tin by lead and sulfur by selenium 
substitution on crystals properties is discussed. 



I. METHOD OF CALCULATIONS 

The calculations of the band structure of both 
phases of Sn2P2S6 crystal as well as energy levels of 
P2S6 molecule has been performed by means of the pack- 
age program ABINIT— (total and projected densities of 
states where calculated using SIESTA 4 ^ software) in the 
framework of the DFT, using the local density approxi- 
mation for representing the exchange-correlation inter- 
action. A basis set of 28000 plane waves, restricted 
by the kinetic energy E cut = 25 Hartree has been used. 
The tin, sulphur and phosphorus atoms had the following 
electron configurations: Sn: [Kr] 5s 2 5p 2 ; S: [Ne] 3s 2 3p 4 ; 
and P: [Ne] 3s 2 3p 3 , respectively. The 'frozen' core elec- 
tron configurations for each atom is shown in brackets. 
The first-principles pseudopotentials in the Hartwigsen- 
Goedecker-Hutter scheme^ 2 - have been applied and the 
integration over irreducible part of the Brillouin zone 
has been done by means of the thetrahedron method us- 
ing the 4x4x4 Monkhorst-Pack mesh^ 3 - of k-points. 
The chosen parameters were sufficient for a good conver- 
gency in the calculations. Prior to commence the calcula- 
tion of physical properties of Sn2P2S6 crystal, we carried 
out the structural optimization, which minimized total 
energy of the system simultaneously with the forces 44 
acting on atoms. The spin-orbit interaction was not 
taken into account in our calculation. The parameters 
of Sn2P2S6 crystal obtained after structural relaxation 
can be compared with experimental data, presented in&2. 
The comparison of experimental and calculated relaxed 
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values of lattice constants demonstrate the difference of 
about 3% in order of magnitude. 



II. ELECTRONIC STRUCTURE OF P 2 S 6 
CLUSTER 

The molecular orbitals of P2S6 cluster create cova- 
lent P — S and P — P bonds. Their hybridization 
with tin atomic orbitals determines electronic structure 
of Sn 2 P2S6 crystal. The electronic energy spectra of 
this material could be analyzed by calculation of free 
P2S6 cluster electronic structure with following account- 
ing of their molecular orbitals hybridization with atomic 
orbitals of tin. The calculated energy spectrum and par- 
tial densities of electron states for s, p and d orbitals of 
phosphorus and sulfur atoms of free P2S6 molecule il- 
lustrate the formation of different molecular orbitals at 
creation of covalent P — S and P — P bonds. The spatial 
electron density distribution for related energy levels re- 
flects peculiarities of these bonds. It is seen that for clus- 
ter the energy level near —15 eV mostly is determined by 
hybridization of phosphorus s orbitals. Here it have been 
also found some contributions of sulfur s and p orbitals. 
The hybridization of these atomic orbitals (scheme (1) at 
(Fig. [2])) creates bonding P — P and P — S molecular 
orbitals. The level near —14 eV is determined by anti- 
bonding combination of two phosphorus s orbitals and 
by bonding hybridization of s orbitals (2) of phosphorus 
and sulfur atoms. 

The group of the levels near —12 eV mostly is formed 
by sulfur s orbitals. These levels are related to the P - 
S bonding and P — P antibonding molecular orbitals 
(3,4,5,6). The levels near —6.5 eV and —8.5 eV appear as 
a replica of doublet of the levels near —15 eV and —14 eV 
in results of antibonding hybridization of phosphorus and 
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FIG. 2. The partial densities of states and space distribution 
of electron density for the molecular orbitals of P2S6 cluster 
in free state. 



sulfur atoms s orbitals. The molecular orbitals for levels 
near —8.5 eV are P — P bonding and P — S antibonding 
(7), the orbitals with energy near —6.5 eV are antibond- 
ing for all P — S and P — P bonds (8). In interval 
between —3.5 and —4.7 eV, the energy levels are created 
by bonding hybridization of phosphorus and sulfur p or- 
bitals. These orbitals (9, 10, 11, 12, 13) are bonding for 
P — P and P — S bonds. In region from till —1.5 eV, 
the energy levels of cluster are also formed by p orbitals 
of phosphorus and sulfur atoms. They are hybridized in 
P — P bonding and P — S antibonding molecular or- 
bitals (14 -23). Here the contribution from phosphorus 
d orbitals is also presented. 

It is worth to be noted that upmost 23rd energy level 
of P2S6 molecule is double degenerated due to high sym- 
metry but only half occupied. Such peculiarity can cause 
instability of this complex by means of Jahn- Teller like 
mechanism. 

From detailed analysis of the electron energy spectrum, 
the important information about creation and character 
of the chemical bonds in P2S6 cluster could be found 
(Fig. [2]). The bond P — P is determined by a hybridiza- 
tion of phosphorus s orbitals (the level near —15 eV (1)) 
and by their replica near —8.5 eV (7). In P — P bond, the 
contribution from tt hybridization of p x and p y orbitals of 
phosphorus, that are oriented normally to the bond direc- 
tion, is also presented. The levels of these orbitals (9) are 
placed near —4 eV. However, the essential contribution 
into energy of the P — P bond adds the a hybridization of 
phosphorus p z orbitals which are oriented along the bond. 
The bonding combination a (p z + p z ) (22) has filled by 
electrons with the energy level in the range -. — 1.5 eV. 
Also the nonbonding combination cr* (p z — p z ) which is 
related to the empty energy level has been found in this 
energy region. Some contribution to the bond between 
PS3 structure pyramids of P2S6 cluster also came from 
hybridization of p orbitals of sulfur atoms that belong to 
different pyramids. 

By a hybridization of phosphorus and sulfur s orbitals, 
the P — S bonds are created. Also the hybridization of 
p orbitals of phosphorus and sulfur atoms is observed 
(Fig. [2]). Such hybridization has obviously both a and tt 
character. Thus the next scheme for appearing of molec- 
ular orbitals that form P — S bonds in PS3 structural 
pyramid could be proposed. The phosphorus atom re- 
alizes sp 2 hybridization from which three symmetrically 
oriented bonds with involving of sulfur p orbitals appear. 
Two electrons from surrounding cations in a crystal lat- 
tice and one s electron of phosphorus atom (excited on d 
orbital) supply filling of covalent P — S bonds. 

In Sn2P2S6 crystal, all energy levels for P2S6 clusters 
are occupied, and in the ionic bonding approach, the 
charge states S 2_ for sulfur ions and P 4+ for phosphorus 
ions are expected. However, as it follows from calcu- 
lations for P2S6 molecular orbitals, enough high charge 
density is found at phosphorus atoms. By this matter, 
not high positive charge (drastically smaller than +4) is 
expected for the phosphorus ions. In addition, the cal- 
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FIG. 3. The hybridization scheme for electronic orbit als in 
Sn 2 P2S6 crystal 



dilations provide evidence about high electronic density 
at the middle of P — P bond. 



III. ELECTRONIC STRUCTURE OF Sn 2 P 2 S 6 
CRYSTAL 

At building of Sn2P2S6 crystal structure with two for- 
mula units in the elementary cell (Fig. [1]), the quantity 
of energy levels of P2S6 clusters is doubled with their 
energies splitting in results of the inter-cluster interac- 
tion. Also, the energy levels of 5s 2 orbitals of four Sn 2+ 
cations are added to the structure of the crystal valence 
band (VB). The cation's bp 2 orbitals participate in for- 
mation of the conduction band of crystal. The scheme of 
electron orbital hybridization in the crystal could be pre- 
sented as additive combination of above described scheme 
of P2S6 cluster orbitals creation with scheme of these 
molecular orbitals hybridization with tin atomic orbitals. 
Generally, for the crystal, four tin atoms, four phospho- 
rus atoms and 12 atoms of sulfur in the elementary cell 
have 108 valence electrons that are placed at 54 energy 
levels in the VB. 

It's known^ that energy positions of the atomic or- 
bitals of phosphorus, sulfur and tin are next: P 3p = 
-8.35 eV, P3s = -17.13 eV; S 3p = -10.28 eV, 
S3s = -20.8 eV; Sn bp = -4 eV, Sn 5s = -11 eV. 
The scheme of hybridization (Fig. [3j) could be proposed 
which is in agreement with calculated energy spectra for 
free P 2 S 6 cluster (Fig. [2]) and for Sn 2 P 2 S 6 crystal (Fig.©. 
At the crystal structure formation, the energy of phos- 
phorus valence orbitals almost doesn't change while the 
bonding energy for sulfur valence orbitals strongly lowers 
(almost by 7 eV). This is in agreement with raised elec- 
tronic density on the sulfur anions and with enough high 
electronic density surrounded the phosphorus atoms. 

In accordance to the calculated energy spectra and 
densities of states of Sn2?2S6 crystal (Fig. 2]), its VB 
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FIG. 4. The electron energy spectrum of Sn2P2S6 crystal in 
the paraelectric (top) and the ferroelectric (bottom) phases. 
The calculated total density of states is compared with ex- 
perimental XPS data££. 



could be divided into eight subbands which are labeled 
at the hybridization scheme of atomic and molecular elec- 
tronic orbitals. Remember, that for the energy spec- 
trum of free P2S6 cluster only seven subbands were found 
(Fig. [2]). For Sn2?2S6 crystal, the additional levels of tin 
atomic s orbitals are placed near — 8 eV and they are 
related to the sixth subband of the VB. 

The subband I contains the two energy levels near 

— 16.5 eV for which a contribution of phosphorus s or- 
bitals is dominated (Fig. [5]). Here, a small contribution 
of sulfur s orbitals is also presented, and a minor appear- 
ance of sulfur p orbitals is observed. The contribution of 
tin atoms valence orbitals in this subband is specific pe- 
culiarity of the crystal energy spectrum. Generally, the 
bonding orbitals of covalent P — P and P — S bonds are 
created in subband I (Fig. [6]). 

The subband II includes the two energy levels near 

— 15.5 eV. It is formed by antibonding combination of two 
phosphorus atoms s orbitals and by bonding hybridiza- 
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FIG. 5. The total and partial electron densities of states for 
Sn2P2S6 crystal in the paraelectric (denoted by p) and ferro- 
electric (denoted by f) phases. 



tion of s orbit als of phosphorus and sulfur atoms. For 
these levels also, a some contribution from phosphorus d 
orbitals is observed. 

The subband III in region from —13 eV to —13.8 eV 
has eight energy levels for which electron charge den- 
sity is mostly localized at sulfur atoms (Fig. [6]). Here 
the contribution of phosphorus 3p and Ad orbitals is also 
presented. This subband is characterized by bonding hy- 
bridization for P — S and antibonding hybridization for 
P — P covalent bonds in the crystal structure. 

The subbands IV and V near —10 eV and —8.5 eV 
(both of them contain two energy levels) are formed by 
s and p orbitals of phosphorus and sulfur atoms. They 
are the replica of the subbands II and I and originated 
from their hybridization with subband III. For the IV-th 
subband, the charge is mostly localized at phosphorus 
atoms and it has bonding character for the P — P bonds 
and antibonding character for the P — S bonds. The 
subband V has antibonding character for both P — P 
and P — S bonds. 

The subband VI with four energy levels in range from 
—8 up to —6.5 eV appears in the VB of Sn2?2S6 crystal 
as a result of hybridization of tin atomic electron orbitals 
with P2S6 clusters molecular orbitals. This subband is 
mainly formed by tin s orbitals and by p orbitals of phos- 
phorus and sulfur (Fig. [5]). Charge of this subband's hy- 
bridized orbitals is located between phosphorus atoms 
and around tin atoms and has Sn — S and P — P bond- 
ing character (Fig. [6]). The peculiarities of orbitals hy- 
bridization for every of four levels from this subband will 
be analyzed in details later with the aim of interatomic 
interactions explanation that are related to the ferroelec- 
tric phase transition in Sn2P2S6 crystal. 

The subband VII located in energy region between 
—6.5 and —4.5 eV and contains ten energy levels. They 
are formed by phosphorus and sulfur p orbitals and also 



include small contribution of tin electronic orbitals. This 
subband has bonding character for P — S and P — P 
bonds and it is antibonding for the Sn — S bonds. 

The VHI-th subband is situated near the top of crys- 
tal's valence band and includes 24 energy levels in the 
energy range from —4.5 till —0.5 eV. The considered sub- 
band is mostly formed by lone pairs of sulfur p orbitals, 
with some participation of phosphorus p orbitals, and it 
has P — P bonding and P — S antibonding character. 
Here the hybridization of s and p orbitals of tin atoms 
what determines their stereoactivity (Fig. [5]) is also re- 
produced. The nature of Sn 2+ cations stereochemical ac- 
tivity in the Sn2P2S6 crystal structure will be analyzed 
in details later on. 



IV. HYBRIDIZATION OF TIN ATOMIC 
ORBITALS WITH MOLECULAR ORBITALS OF 
P 2 S 6 CLUSTERS 

As it was mentioned above, the valence subband VI 
includes four levels in the region from —8 to —6.5 eV 
which are related to s orbitals of tin atoms (Fig. HH6]). 
These are levels from 17 to 20, and their spatial electronic 
charge distribution illustrates a bonding peculiarities of 
tin atomic orbitals with sulfur and phosphorus orbitals 
which create the clusters P2S6. The level 17 is charac- 
terized by enough strong bonding of tin atoms with two 
nearest sulfur atoms, the levels 18 and 19 demonstrate 
strong bonding of tin atom with one of the nearest sulfur 
atom. The spatial charge distribution of the 20th level 
is of special interest - here the electron density is elon- 
gated from tin atom to the middle of P — P bond inside 
P 2 S 6 cluster (Fig. [7^). 

It can be easily seen that in the paraelectric phase in 
addition to — Sn — S — P — P — S — Sn — sequence 
of chemical bonds the — Sn — P — P — Sn — series also 
exists. Such sequence of direct bonds of tin atoms with 
phosphorus atoms appears due to anisotropy of the spa- 
tial charge distribution for the level 20 (Fig. 03). Named 
distribution has a form of layers that are oriented close to 
the plane (101). It is important that directions of the tin 
atoms shifts at transition into ferroelectric phased are 
also rather close to orientation of the mentioned plane. 

In the ferroelectric phase, the two pairs of nonequiva- 
lent tin atoms are presented in the crystal structure what 
is clearly illustrated by the spatial distribution of the 
electron charge density for the 20th level (Fig. [8^). Near 
mentioned P2S6 cluster, one of the tin atoms, further 
Sn2, is approached to the middle of P — P bond. At 
this, the electron charge distribution which surrounded 
the Sn2 atom connects with charge distribution between 
two phosphorus atoms. Another tin atom (Sni) goes 
away from the middle of P — P bond, and surrounded 
Sni atom electron charge distribution is oriented to one 
of the sulfur atoms. Such difference in hybridization of 
the electron orbitals for two types of tin atoms in the 
ferroelectric phase determines disappearance of the layer- 




FIG. 6. The spatial charge density distribution (in plane which contains S 
subbands I - VIII for the paraelectric phase of Sn2P2S6 crystal. 



S bonds) for orbitals in the valence 




FIG. 7. The electron density distribution - (a) - in plane 
parallel to P — P bond, (b) - in perpendicular one - for 
orbitals of 20th energy level in the valence subband VI for 
the paraelectric phase of Sn2P2S6 crystal. 



like anisotropy for the electron spatial charge distribution 
(Fig. Eh). 



V. TRANSFORMATION OF ELECTRON 
ENERGY SPECTRA AT TRANSITION FROM 
PARAELECTRIC PHASE INTO 
FERROELECTRIC PHASE 

According to experimental datai£~— the chemical 
bonds and electron energy spectra of Sn2P2S6 crystal 
have noticeable changes at the ferroelectric phase tran- 
sition. The calculated energy spectra demonstrates the 
changes in the energy gap and in positions of all energy 
levels of the VB - energies of electron states density peaks 




FIG. 8. The electron charge distribution - (a) - in plane paral- 
lel to P — P bond, (b) - in perpendicular one - that illustrate 
the hybridization of tin and phosphorus orbitals for the 20th 
energy level in the valence subband VI of the Sn2P2$6 ferro- 
electric phase. 



in the VB shifts approximately by 0.5 eV (Fig. I4|5|) . In 
acentric phase the degeneration of electron energy levels 
disappears in some regions of the Brillouin zone, what 
determine higher smearing of the energy distribution of 
electron density of states. It is important to mention that 
the lowering of the electron density of states near the top 
of the VB is also observed at transition into ferroelectric 
phase together with a rise of the energy gap. 

The changes of electron energy spectra obviously re- 
flect an important role of electron-phonon interaction in 
nature of the Sn2P2S6 crystal spontaneous polarization. 
Such interaction is illustrated by a transformation of 
the spatial electron density distribution at change of the 
atoms coordinates in the crystal structure. The squares 
of wave functions for electron orbitals, which are summed 
for the energy levels of valence subband VI in the para- 
electric and the ferroelectric phases of Sn2P2S6 crystal, 
are shown at Fig. [9] In the elementary cell of centrosym- 
metric structure, a similar distribution of electron den- 
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sity around four tin atoms is observed. This distribution 
reflects the stereoactivity of the lone electron pair 5s 2 
of cations Sn 2+ . Also, the spatial charge distribution 
around phosphorus atoms is similar. This is in agree- 
ment with the presence of inversion center at the middle 
of P — P bond. In the acentric structure, the two pairs 
of tin atoms, with different distribution of surrounding 
charge, appear. The nonequivalence of the electron den- 
sity distribution near the phosphorus atoms is also seen. 
It should be mentioned that elevated electron density is 
located in the vicinity of neighboring tin atoms Sn2 and 
phosphorus atoms Pi. Also, lowering of the surround- 
ing charge is observed for the adjacent atoms Sni and 
P2. Acentricity of the P2S6 clusters is also reflected in 
deformation of the charge distribution along the P — P 
bonds. 

Let's analyze peculiarities of the orbit als hybridization 
for study of the tin 5s 2 electrons stereoactivity and their 
contribution into lattice spontaneous polarization. In the 
crystal structure the tin atoms are placed in the polyhe- 
drons that are created by eight sulfur atoms. At be- 
ginning we will consider the mixing of tin orbitals with 
valence orbitals of surrounding sulfur atoms. Further, 
it will be analyzed a role of hybridization between tin 
atomic orbitals and P2S6 cluster molecular orbitals in 
change of relation between short range and long range 
interactions in Sn2?2S6 crystal lattice, which induce fer- 
roelectric phase transition. 



VI. STEREOACTIVIY OF LONE 5s 2 
ELECTRON PAIR OF Sn 2+ CATIONS 

By first principles calculations in approach of frozen 
phonons, it was founcU^ that in Sn2P2S6 crystal the polar 
optic mode with symmetry B u in the paraelectric phase 
could be destabilized only at account of their nonlinear 
interaction with fully symmetrical A g optic mode. Gen- 
erally, it is necessary to include the nonlinear interaction 
of A g B 2 type for all 13 B u and 15 A g normal vibrations 
of the lattice with symmetry P2\/c in the paraelectric 
phase. Thereafter it could be then explained the varia- 
tion of the atomic coordinates at transition into the ferro- 
electric phase with symmetry Pc - two tin cations (that 
are related by the symmetry plane) have some shift rela- 
tively to the anion sublattice, and another two tin cations 
found a bigger shift (flipping) relatively to their positions 
in the paraelectric phase. It was determined, that effec- 
tive potential in the normal A g -B u coordinates has three 
minima - central one reflects the metastable paraelectric 
state, and two side minima are related to two domains of 
the ferroelectric phase in uniaxial ferroelectrics. 

The complicate potential relief for Sn2P2S6 crystal is 
obviously determined by relaxation of electron lone pair 
of tin cations. The electron pair in 5s 2 configuration is 
stereoactive - s orbitals of cations are hybridized with 
p orbitals of sulfur. The antibonding energy level as 
a result of such hybridization is still occupied by elec- 



trons. For the energy gain, this level interacts with tin 
p orbitals - sp 2 hybridization is realized for which the 
bonding level lowers their energy and antibonding level 
elevates higher of Fermi level. Such hybridization is pro- 
portional to the acentricity of surrounding crystal field, 
and it induces deformation of built by eight sulfur atoms 
dodecahedron which surround the tin cation. Thereafter 
the metal cation shifts away from the dodecahedron cen- 
ter, and space distribution of the electron charge is clearly 
different from spherical. The stereoactivity of electron 
lone pair represents the second order Jahn- Teller effect - 
SOJT. 

Already in the paraelectric phase, the tin cations have 
placed in general position of the elementary cell and their 
surrounding by electron density is definitely not spheri- 
cal - almost rigid dipoles exist which are not aligned. At 
cooling the stereoactivity of electron lone pair growths 
and the orientation of cation shifts in nearest elementary 
cells is correlated by dipole-dipole interaction. The spon- 
taneous polarization appears with two contributions - 
"displacive" and " order /disorder" . 

The experimental data of XPS spectroscop y 20 ! 30 about 
electron energy spectra near the top of VB in the para- 
electric and the ferroelectric phases of Sn2P2S6 crystal 
confirm the lowering of the electron energy states density 
at transition into acentric structure. The calculations of 
energy spectra show big enough contribution of tin s or- 
bitals into electron states density near the top of the VB 
and permit a possibility for tracing of their change at 
the phase transition (Fig. [5j). The calculated space dis- 
tribution of total electron density illustrates an evolution 
of electron lone pair at transition from the paraelectric 
phase into the ferroelectric one (Fig. [TP]) . The appear- 
ance of the tin cations nonequivalence is observed: two 
of them shifts in a direction of those sulfur atoms with 
which in the paraelectric phase they have the biggest 
overlap of electron orbitals; alternative two cations "flip" 
the biggest overlap of their electron orbitals in direction 
of other sulfur atoms (Fig [TO]). 

The high coordination of tin cations could be obviously 
described by taking into account their d orbitals. For 
symmetry requirements the hybridization of sp 2 d 5 type 
satisfies and they could be related to the positioning of 
the tin cations inside of the eight caped polyhedron of 
sulfur anions. Indeed, the orbitals s, p Xl p Zl d XZl d xyi 
d Z y, d z 2, d z 2_ x 2 transform on the irreducible presenta- 
tions A g , B U1 B U1 A g , B g , B g , Ag, Ag (orientation of 
the monoclinic symmetry plane coincides with cartesian 
plane XZ) that have fully symmetrical combination. 

In the elementary cell of the paraelectric phase, the 
four such dodecahedrons are equivalent - they pair-by- 
pair are related by second order screw axis or by glade 
mirror plane (Fig. [TTa). In the ferroelectric ph ase as a 
result of earlier discussed charge density redistribution, 
the inversion center and symmetry axis disappear. The 
two pairs of nonequivalent dodecahedrons appear. In one 
type of the dodecahedrons at the phase transition, the 
electron density switches between two nearest sulphur 



FIG. 9. The electron density space distribution for the VI valence subband in the paraelectric (a) and the ferroelectric ( 
phases of Sri2P2S6 crystal. 



atoms which corresponds to the strongest Sn — S bond. 
In the other type of dodecahedrons at cooling from the 
center symmetric phase to acentric one, the "flipping" of 
the electron density between almost oppositely oriented 
Sn — S bonds is observed (Fig. HTb). 

The built by sulfur atoms pyramids could be divided 
in the coordination polyhedrons (Fig. [TT]) . Localized at 
hybridized sp 2 -\ike orbitals of tin atoms, the electron den- 
sity is oriented in direction opposite to the base of named 
pyramids. This is a direction to the three sulfur atoms 
with weaker Sn — S bonds. 

A characterization of the chemical bonds changes at 
phase transition could be found by compare of calcu- 
lated data (Mulliken charges and overleap parameters 
for the electron orbitals of neighbor atoms) with ex- 
perimental structure data, with Mossbauer, XPS and 
NMR spectroscopies data. For the paraelectric phase 
of Sn2P2S6 crystal, the following electron configura- 
tions were calculated: for four equivalence tin atoms 
Sn - 5s 1,865 5p 1,153 5<i 0,222 ; for four equivalence phospho- 
rus atoms P - 3s 1 187 3p 2 - 508 3d 1172 ; for one of the three 
types of sulfur atoms S - 3s 1833 3p 4 - 231 3d 0179 . It could 
be seen that d orbitals are populated, what agrees with 
the explanation of high dodecahedral coordination of tin 
cations on background of sp 2 d 5 type hybridization. The 
highest population of d orbitals is appeared for phospho- 
rus atoms. As was mentioned earlier, for the P2S6 clus- 
ters the bonds in PS3 structure pyramids are deter- 
mined by sp 2 hybridization of phosphorus atomic orbitals 
with their further a hybridization with sulfur p orbitals. 
Thereafter one of two 3s 2 electrons of phosphorus is in- 
volved into sp 2 hybridization of atomic orbitals, other 
electron is excited on the atomic d orbital. The named 
sp 2 hybridization involves two p electrons of phosphorus, 
third p electron by <j(p — p) hybridization creates cova- 
lence P — P bond. 

In the ferroelectric phase, the electron configurations 
for two types of tin atoms are the following: Sni - 



5s 1.851 5p 1.163 M 0.224. _ ^1.8505^1.177^0.226. It ig 

seen that at transition into the ferroelectric phase, the 
quantity of s electrons lowers by An s = 0.014. Such cal- 
culations data are in agreement with observed lowering of 
isomer shift for spectral line of 119 Sn Mossbauer effect at 
cooling from the paraelectric phase into the ferroelectric 
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The diminishing of s electron quantity in the ferro- 
electric phase for two types of tin atoms is compensated 
in a different step (level) by growing of p orbitals oc- 
cupancy. Hence the total charge of tin atoms changes 
from 3.240e in the paraelectric phase to 3.239e (Sni) and 
3.251e (Sn2) in the ferroelectric phase. The calculated 
rise of electron density in the vicinity of Sn2 atoms is 
in agreement with decreasing of resonance frequency in 
NMR spectrum for 119 Sn from —781.3 ppm in the para- 
electric phase to —782.5 ppm in the ferroelectric phased. 
The increase of the resonance frequency till —754.5 ppm 
for Sni is obviously determined by a occupancy redistri- 
bution from s orbitals to p orbitals at some lowering of 
total charge. 

Determined by positions of XPS spectral lines 2 ^, the 
energy of chemical binding for the Ad core orbitals in the 
ferroelectric phase differently increases for the two types 
of tin cations. This is obviously defined by growing of 5s 2 
electrons lone pair stereoactivity and by these electrons 
redistribution on more distanced from tin cores p orbitals. 

The changes in Sn2?2S6 crystal structure at the spon- 
taneous polarization appearance are characterized by the 
calculated values for the overleap parameters of electron 
orbitals of tin atoms and the nearest sulfur atoms. In 
the paraelectric phase such parameter with the biggest 
value about 0.054 has been found for the shortest bonds 
Sn — S in the dodecahedrons of sulfur atoms (Fig. [TT]) . 
It must be mentioned that in the nearest dodecahedrons, 
which are related by a second order screw symmetry 
axis, such Sn — S bonds have opposite orientations of 
their projections onto plane (010), which contains the 
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FIG. 10. The transformation of the electron density space 
distribution in the vicinity of tin cations at the ferroelectric 
phase transition in Sn2P2S6 crystal. The strongest Sn — S 
bonds is shown by dashed lines. 



spontaneous polarization. In the ferroelectric phase the 
nonequivalence of two pairs of cations Sni and Sn2 and 
surrounded them dodecahedrons of sulfur ions appear. 
Here the biggest overleap parameter for atomic orbitals 
(0.096) was found for the Sn2 — S bond. This bond 
was strong already in the paraelectric phase and it is 
in addition enforced at spontaneous polarization appear- 
ance - the atom Sn2 in the ferroelectric phase is shifted 
in the direction of sulfur atoms with the biggest content 
of chemical bonds covalency. For other coordination do- 
decahedron the strongest chemical bond Sni — S has 
the overlap parameter 0.087. Here at transition from the 




FIG. 11. The transformation of sulfur atoms coordination 
polyhedrons around the tin cations at the ferroelectric phase 
transition in Sn2P2S6 crystal. The sulfur atoms that have 
stronger bonds with tin atoms are linked by dark lines. The 
strongest Sn — S bonds is shown by dashed lines. 



paraelectric phase into the ferroelectric one, the "flip- 
ping" of maxima in the space distribution of electron 
density occurs between bonds of central metal atom and 
ligand atoms in the coordination polyhedron. As a result 
of such "flipping" in both types of the dodecahedron, the 
strongest bonds Sni — S and Sn 2 — S have nearly ori- 
ented projections onto plane (010) fFig. [TTb). 

The changes in the electron density space distribution 
correlate with variations of the interatomic distances. For 
example, in the paraelectric phase the Sn — S bonds 
with the biggest overleap parameter (0.054) are strongest 
and they have the smallest length in the ferroelectric 
phase. For the dodecahedrons with tin atoms of Sn 2 type, 
the overleap parameter for strongest bonds increases to 
0.096 and their length decreases by 0.21 A. Hence the 
occupancy of XZ plane oriented p orbitals raises - from 
the calculations follows: Ap x = 0.048, Ap y = —0.002, 
Ap z = 0.022. For the dodecahedrons with tin atoms 
of Sni type at the electron charge density "flipping" 
on almost oppositely oriented Sn — S bond, the over- 
leap parameter, for the strongest bond in the paraelec- 
tric phase, lowers from 0.054 till 0.022 and length of 
this bond growths by 0.26 A. The strongest Sni — S 
bond in the ferroelectric phase is characterized by over- 
leap parameter 0.087 and its length decreases by 0.35 A. 
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Here the changes for occupancies of tin p orbit als are: 
Ap x = 0.025, Ap y = 0.036, Ap z = 0.010. We could 
see that cations of Sn2 type donate the biggest electronic 
contribution into the spontaneous polarization and this 
contribution has "displacive" character. For the cations 
of Sni type, the electronic contribution is a little smaller 
and this one has "ordering" character. 

The calculated changes of Mulliken charges and atomic 
orbitals overleap parameters coincide with the temper- 
ature dependence of the resonance frequencies in NMR 
spectra for 31 P phosphoru a 18 i 19 . In the paraelectric phase 
all phosphorus atoms are equivalent (calculated charge 
is 4.869e) and here only one NMR line with frequency 
92.12 ppm is observed. At the spontaneous polarization 
appearance, the inversion center vanishes what is asso- 
ciated with growth of phosphorus atoms nonequivalence 
in P2S6 clusters. For two types of phosphorus atoms the 
following electron configurations were calculated: Pi - 

35 1.1909 3p 2.499 3d 1.183 ; p 2 _ 3^1.1823^2.5093^.193^ For 

atoms of Pi type with grown calculated charge (3.251e) 
in result of shielding effect, the NMR specter resonance 
frequency decreases to 89.2 ppm. This is in agreement 
with closeness of Pi type atoms and Sn2 cations for which 
also a lowering of the NMR resonance frequency is ob- 
served^ in the result of electron density growing in their 
vicinity at cooling into the ferroelectric phase. For atoms 
of P2 type in the ferroelectric phase, the resonance fre- 
quency rises till 93.7 ppm. In this case the weakening 
of the charge shielding effect is obviously determined by 
diminishing of the s orbitals occupancy. By the way, 
for neighbor Sni cation the quantity of s electrons also 
decreases, what induces growth of related resonance fre- 
quency in 119 Sn NMR spectrum^. 



VII. DISCUSSION OF RESULTS 

The calculated electron energy spectra, densities of 
electron states and their variation at transition from the 
paraelectric phase into the ferroelectric one coincide with 
the available structure data and results of experimen- 
tal investigation of the chemical bonds nature. On this 
ground the next generalized description of sources of the 
spontaneous polarization appearance in Sn2P2S6 crystal 
could be proposed. 

The ferroelectric distortion is proportional to the dif- 
ference between the Sni and Sn2 positions. For this dis- 
tortion the short range (mostly Sn — S) repulsions must 
be sufficiently small in order to allow to shift the equi- 
librium Sn positions from the center of the chalcogen 
dodecahedron. The effective charge of phosphorus has to 
be sufficiently small in order to allow the shift of the Sn 
cations in the direction of P — P bond. Such require- 
ment could be satisfied by the following way. At fully 
symmetrical A g lattice vibration, the important changes 
of electron charge distribution in the elementary cell oc- 
cur - the charge partially is waded from anions P2S6 onto 
cations Sn. Hence the p orbitals of tin cations have to be 



occupied - the stereoactivity of valence electrons of these 
cations is realized by their partial hybridization with p or- 
bitals of neighbor sulfur atoms. Such hybridization low- 
ers the short range repulsion between tin cations and the 
nearest sulfur atoms what govern their approaching. The 
charge transferring between tin and sulfur atoms man- 
ages some lowering of electrostatic interactions energy. 

Generally, in the ground state (at K) the metastable 
center-symmetric Sn2P2S6 structure is possible, for 
which both opposite tin atoms are a little approached 
to the middle of the P — P bond. However, at low 
temperatures the acentric structure is energetically more 
favorable. At approaching of two tin cations, that are 
related by inversion center at middle of P — P bond, 
their Coulomb repulsion increases (such repulsion be- 
tween cations of tin and phosphorus obviously don't play 
important role because positive effective charge of phos- 
phorus is not big). Hence energetically more advanta- 
geous could be approach of one tin cation to the center 
of P2S6 cluster at the repulsion of opposite tin cation. In- 
deed, in the center-symmetric structure both tin atoms 
are placed on distance 3.633 A far from middle of P - 
P bond, in acentric structure such distances equal to 
3.463 A for Sn2 atom and 3.857 A for Sni atom. 
The space between considered tin atoms increases from 
7.266 A in the paraelectric phase till 7.310 A in the 
ferroelectric phased. 

Removing of Sni atom away from P2S6 cluster de- 
creases hybridization of his valence electron orbitals with 
molecular orbitals of cluster. By this matter the electron 
charge in cluster moves onto PS3 structural pyramid with 
Pi atom at their top that is the nearest to Sn2 atom. 
At repulsion of Sni atom from the cluster, an important 
change of the hybridization character occurs, which is ac- 
companied by localization of valence electrons near the 
tin ion core and by growth of their kinetic energy. Such 
processes have obviously activation character and they 
determine presence of the energy barrier between central 
and side minima in the three-well potential. 

By such way, the induced by the fully symmetrical A g 
vibration, the important changes of charge gradient in 
the elementary cell determine variation of the electron 
configuration for the ions of crystal lattice. The recon- 
struction of electron configuration modifies the balance of 
interatomic interactions what induce instability of B u po- 
lar lattice vibration. Exactly by such manner, the mech- 
anism of lattice modes A g B\ nonlinear interaction could 
be presented, and this one governs the tree-well poten- 
tial presence for fluctuations of the order parameter of 
the ferroelectric phase transition in Sn2P2S6 crystal. 

The value of energy barrier in the three-well poten- 
tial^ equals near 0.015 eV, and the energy difference 
between central and side minima is near 0.01 eV. Such 
energetic characteristics are in agreement with our calcu- 
lation of the electron energy spectra of Sn2P2S6 crystal 
in the paraelectric and the ferroelectric phases. Thus, 
for the ferroelectric phase the full energy was found by 
0.0078 eV smaller in compare with calculated full energy 
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for the paraelectric phase. It must be mentioned that at 
the crystal symmetry lowering the negative contribution 
of coulomb interactions into full energy decreases. 

The appearance of spontaneous polarization in ferro- 
electric crystal is determined by variation of the chemical 
bonds covalency and by delicate balance between short 
range repulsion forces, which determine the relief of local 
potential for the phase transition order parameter, and 
long range displace forces that define energy of intercell 
interaction. 

The hybridization of tin and sulfur atomic orbitals 
which defines the appearance of "partially rigid" electric 
dipoles (pseudospins) in result of tin cations valence elec- 
trons stereoactivity was early analyzed. Such hybridiza- 
tion could be described as sp 2 or sp 2 d 5 combination of 
tin and sulfur atomic orbitals. This fact is clearly demon- 
strated by presence of enough high density of Sn 5s states 
near the top of VB. In addition to that, it have been 
found neediness of accounting for hybridization between 
tin atomic orbitals and P2S6 molecular orbitals what has 
argument by presence of 5s and 5p states of tin even at 
the bottom of VB - in the energy range near —17.5 eV 
where the s states of phosphorus atoms are dominated. 

The obtained pictures of the electron density space dis- 
tribution show (Fig.[7j) the presence of — Sn — S — P — 
P — S — Sn — short-range bonds chains in the paraelec- 
tric phase. The occurrence of the short-range interactions 
in chains of — Sn — P — P — Sn — type also have been 
demonstrated. Such sequences of the short-range inter- 
actions together with long-distant coulomb interactions 
determine the mean-field, which induce correlation of the 
pseudo-spins at lowering of disordering influence of heat 
energy. At temperature down till K, the chains of cen- 
tral symmetric structure groups — Sn — P2S6 — Sn — 
, which are related to the pseudo-spins position in the 
central well of local potential, could obviously exist also. 
However, the correlated ordering of the structure motives 
("dimers") like — Sn — P2S6 — , which responses to the 
pseudo-spins standing in one of the side well of local po- 
tential, is energetically more favorable. 

It is important to remark about rise of anisotropy of tin 
and phosphorus atoms surround in the ferroelectric phase 
(Fig. [8]) what support increasing of their dynamic or Born 
effective charges. The dynamical transfers of charge are 
expected to be larger when such a hybridization involves 
d states, for which the interactions parameters with other 
orbitals are particularly sensitive to the interatomic dis- 
tance^. Also, the amplitude of Born effective charges is 
not monitored by a particular interatomic distance but 
is dependent on the anisotropy of the Sn environment 
along the — Sn — S — P — P — S — Sn — chains. In 
the paraelectric phase, the S 3p electrons are obviously 
widely delocalized and dynamical transfer of charge can 
propagate along the — Sn — S — P — P — S — Sn — 
chains. In the ferroelectric phase, these chains behave as 
a sequence of "dimers" — Sn — S — P — P — S — 
... for which the electrons are less polarizable. 

The anomalously large dynamical charges produce big 



LO-TO splitting for the ferroelectric soft phonon mode 4 ^. 
Moreover, this feature is associated with the existence of 
an anomalously large destabilizing dipole-dipole interac- 
tion, sufficient to compensate the stabilizing short range 
forces and induce the ferroelectric instability. In mate- 
rials where polar soft modes play a major role, the dy- 
namical charge relate the electronic and structural prop- 
erties^. However, for the Sn2?2S6 crystal the big LO- 
TO splitting for polar modes wasn't observed in phonon 
spectr a 49 ! 50 . Here, at 4.2 K such splitting is in the range 
of 2 -T- 7 cm -1 . At heating to the temperature of phase 
transition in the ferroelectric phase, the LO-TO splitting 
for the lowest energy optic mode of B u symmetry (soft 
mode) reaches only the value of 10 cm -1 . 

The low frequency dielectric susceptibility tempera- 
ture anomaly in Sn2P2S6 crystal don't described only 
by dielectric contribution of the polar lattice vibrations. 
On the data of dielectric spectroscopy 51 in the para- 
electric phase, the dielectric contribution from polar lat- 
tice vabrations into static dielectric susceptibility reaches 
only near ten percent. Obviously in the range of phonon 
frequencies, a significant destabilizing dipole-dipole in- 
teraction doesn't appear. The essential contribution into 
dielectric anomaly appears at frequency lowering into 
submillimeter diapason - here the relaxational dispersion 
have been observed 5 ^, which is obviously determined by 
nonlinear dynamic excitations. 

The above attention was accented on the mixing be- 
tween states of tin valence electrons and orbitals of 
phosphorus and sulfur atoms across all energy range of 
Sn2?2S6 crystal VB. In addition, the enough large den- 
sity of phosphorus s and p states is presented near the 
top of VB also. Such phosphorus orbitals also create the 
conductivity band. The defined facts give evidence about 
enough strong mixing of diffusive s orbitals in structure 
of Sn2?2S6 crystal. Obviously, the effective occupation 
of phosphorus d orbitals, which commonly are enough lo- 
calized, also give evidence of their important role in the 
mechanism of electron-phonon interaction for this ferro- 
electrics. 

For the P2S6 structure cluster, with symmetry D^d in 
free state, the 3s electron orbitals of phosphorus atoms, 
which are placed at bottom of the VB (near —15 eV) and 
their bonding a(p z — p z ) orbitals, which create P — P 
bonds and have energy level near the top of VB (in range 
— 1.5-^0 eV), satisfy the transformation according the A g 
irreducible presentation. In addition to these orbitals of 
free anion cluster, the orbitals with energies near —8.5 
and —4 eV are also involved at formation of P — P bond 
(Fig. El. 

In the Sn2P2S6 structure the symmetry of P2S6 anion 
clusters is lowered, however possibility of effective mixing 
for the wave functions of s and p orbitals of phosphorus 
atoms, which have identical symmetry, is obvious. To 
mentioned series of four sets of P — P bonding orbitals, 
which create high electron density at middle of this bond 
(Fig. [6]), the combination of tin atomic orbital related to 
the energy level 20 in VI subband (in range from —8 till 



12 



—6.5 eV) is added. Obviously, in result of such hybridiza- 
tion the similar changes at the top and bottom of the VB 
are observed at appearance of spontaneous polarization 
in the crystal (Fig. [5j). 

Generally, the high effective charge and large polariz- 
ability of P2S6 anionic clusters together with stereoac- 
tivity of tin electron lone pair determine large electronic 
contribution from all atoms into spontaneous polariza- 
tion of Sn2P2S6 crystal. Such situation is in agreement 
with earlier founded^ involving of all 13 B u modes and 
15 A g modes to dynamic instability of investigated ferro- 
electrics. However, obviously important role also belongs 
to the nonlinear interaction with participation of A u and 
B g nonsymmetric modes. For the point group 2/ra in 
addition to the A g B 2 fully symmetric combination, the 
invariants of A u B g B u type are also present. A signifi- 
cant role of such nonlinear mixing of lattice vibrations 
is reflected in strong internal deformations of P2S6 clus- 
ters at the phase transition. In the crystal electron struc- 
ture, such invariants obviously replicate the hybridization 
of molecular orbitals of P2S6 clusters with participation 
of atomic d orbitals. Such hybridization is clearly illus- 
trated by the spatial distribution of electron density for 
the energy level number 20, which is aligned from tin 
atoms to middle of P — P bond (Fig. [7j). 

The energy decreasing for the s orbitals of phosphorus 
and sulfur at the VB bottom (their contribution is dom- 
inated in the lowest subbands - from I to V), and also 
lowering of energy for tin s orbitals, with their contribu- 
tion across whole range of the VB (Fig. [5]), are obviously 
essential for energetic motivation of the transition into 
the ferroelectric phase. Also some lowering of the p and 
d orbitals energy occurs. 

The electronic structure XPS measurements for 
Sn2P2S6 crystal s 20 ! 30 revealed the chemical shifts of Sn 
and P electronic core states to a higher binding energy 
and of S states to a lower binding energy at the crystal 
lattice formation. This shift suggests a charge transfer 
from Sn and P to S atoms. The binding energies and 
chemical shifts strongly change at the phase transition. 
In the ferroelectric phase, the chemical shifts of Sn and P 
atoms are higher while for S atoms they are smaller. So, 
for all atoms of crystal structure at transition into the 
ferroelectric phase the binding energy for core orbitals 
increases. These data give evidence about localization 
of electron charge in space between atoms, or about en- 
hancement of chemical bonds covalency. Such variation 
of the core orbitals energy agrees with experimentally 
observed and founded at calculations (Fig. I4|5p transfor- 
mation of the VB structure and they support an energetic 
stability of the ferroelectric phase. 

Since covalency increases, there might be a possibility 
for the drastic collapse of the sulfur ionic size (which is 
related to charge transfer). For the S 2_ ion the ionic ra- 
dius is 1.84 A, the covalent radius ~ 1.02 A 53 . If the size 
of sulfurs were small compared to the allowed space then 
the sulfur atoms would be weakly bounded in the lattice. 
In this case an imbalance between the decreased (due 



to small S radius) repulsive forces and the polarization 
forces, what tend to displace the ion from its position, 
also support structure rearrangement. 

In whole, for Sn2P2S6 crystal at transition from the 
paraelectric phase to the ferroelectric one complex evo- 
lution of electron and phonon spectra occurs, which could 
be presented as sequence of five steps. Evidently as first 
factor assists change of the electron density charge distri- 
bution in the elementary cell by fully symmetric breath- 
ing modes A g . Such redistribution of the electron density 
prompts the stereochemical activity of tin cations elec- 
tron lone pair and produces the covalence bonds of tin 
atoms with sulfur atoms (hybridization of sp 2 d 5 type), 
and also with phosphorus atoms, what could be consid- 
ered as second part. 

As third stage could be considered the weakening of the 
short-range repulsion between cations of tin and phos- 
phorus, in result of their charges lowering, and at sig- 
nificant coulomb repulsion of the nearest (related by the 
inversion center) tin cations. Mentioned the second and 
the third factors represent the nonlinear interaction of 
A g B 2 type, they govern an anisotropy of polar shifting 
of atoms in the elementary cell and define appearance of 
the dipole structure motives ( — Sn — P2S6 — ), which 
are related to the polar normal coordinates of B u sym- 
metry. The fourth important factor is the dipole-dipole 
interaction which correlates orientation of local dipoles 
(pseudospins) and defines appearance of the spontaneous 
polarization in the crystal structure. And, at finish, as 
fifth circumstance must be accounted that all low sym- 
metry modes participate in result of permitted nonlinear 
A u B g B u relation in the structure transformation. Such 
combination of structure deformations mirrors participa- 
tion of phosphorus and sulfur d orbitals in the covalent 
bonds of Sn2P2S6 crystal. 

It is interesting to compare the peculiarities of chemi- 
cal bonds in Sn2P2S6 sulfide and Sn2P2Se6 selenide com- 
pounds, and also in the lead contained Pb2?2S6 crys- 
tal. At first, we will consider the binary compounds MX, 
where M - metal Ge, Sn, Pb, X - chalcogen O, S, Se, Te. 
The stereoactivity of electron lone pair for metal atoms is 
determined by sp 2 hybridization of their s and p orbitals 
with p orbitals of chalcogen atoms. Such hybridization is 
determined by the positions of the energy levels of elec- 
tron states and by width of the related energy bands in 
the crystal structur e) 47 ! 54 . The smallest energy difference 
is present between positions of energy levels of germa- 
nium s orbitals and oxygen p orbitals. Consequently for 
the compound GeO, the largest stereoactivity of 4s 2 elec- 
tron orbitals of Ge is observed 5 -^. At transition from Ge 
to Sn and than to Pb, the energy of chemical binding for 
the s orbitals increases. Thus with transition from O to 
S, and further to Se and Te, the energy of chemical bind- 
ing for their p orbitals decreases. It is expected that the 
hybridization of Ge 4s orbitals and O 2p orbitals is the 
strongest, and hybridization for Pb 6s orbitals and Te 5p 
orbitals is the most weak. However the hybridization is 
also influenced by width of related energy bands in the 



13 



crystal structure. The increase of width for the s and 
p electron states energy bands and their overleap could 
partially compensate increase of the energy distance be- 
tween related energy levels, what produce some level of 
the stereoactivity and covalency of M — X bonds. 

Thus, at transition from Sn2P2S6 to PD2P2S6 the bind- 
ing energy for Pb 6s level increases, what make weaker 
the stereoactivity of the 6s 2 electron lone pair in the do- 
decahedron of sulfur atoms. Obviously by this matter, 
the observed 55 ionicity for the Pb — S bonds is higher in 
compare with Sn — S bonds ionicity. The melting tem- 
perature and the energy gap both rise at substitution Sn 
by Pb. The paraelectric phase in Pb2?2S6 is stable at 
cooling till 4.2 K 7 -. 

The largest stereoactivity of Ge 4s 2 electron lone pair 
in surround of sulfur atoms gives natural explanation of 
absence of Ge2P2S6 crystal structure. The Ge atoms 
couldn't be placed in positions with high coordination of 
sulfur atoms. Obviously introducing into Sn2?2S6 crystal 
impurity of germanium in the charge state Ge 2+ strongly 
elevates the temperature of ferroelectric phase transition 
what was observed by dielectric investigations 24 . Cer- 
tainly the impurity in charge state Ge 4+ will not be 
stereoactive and will not support rise of temperature in- 
terval for the ferroelectric phase existence. 

According to just described tendencies, at transition 
from Sn2?2S6 to Sn2P2Se6 the stereoactivity of Sn 5s 2 
electron lone pair in the dodecahedron of selenium will be 
smaller than in the case of sulfide compound. However, 
for the ternary compounds the ion-covalence bonds Sn — 
S and Sn — Se are modified depending on peculiarities 
of P — S and P — Se bonds in P2S(Se)6 anion clusters. 
On the Mossbauer spectroscopy data 55 at substitution 
of sulfur by selenium, the isomer shift for 119 Sn nucleus 
decreases what directly show on higher covalency of Sn — 
Se bonds. The NMR spectroscopy for 119 S n 19 i 56 shows 
increase of the resonance frequency from —781.3 ppm in 
Sn 2 P 2 S 6 to -608 ppm in Sn 2 P 2 Se 6 . The NMR spectral 
line for 31 P decreases its frequency from 92.12 ppm in 
Sn 2 P 2 S 6 -^ to 28.7 ppm in Sn 2 P 2 Se 6 - 5 ^. These data pro- 
vide evidence about lowering of the electron density in 
vicinity of tin nucleus and about their rise at phosphorus 
nucleus at transition from sulfide to selenide compound. 
The mentioned tendency could be explained by smaller 
electronegativity of selenium. Obviously the bonds P - 
Se are less polar what improve higher electron charge 
of phosphorus atoms in anionic clusters. Such situation 
probably also support more effective hybridization be- 
tween cluster molecular orbitals and tin atomic orbitals 
(mostly the bonding hybridization of tin orbitals with 
molecular orbitals that are localized in middle of P — P 
bond) and increases the stereoactivity of tin lone pair of 
electrons in Sn2P2Se6 crystal. 

The growth of covalency and weakening of electro- 
static interactions determine lowering of the melting 



temperature for Sn2P2Se6 crystal in compare with 
sulfide analog, define decrease of the energy gap and 
govern smaller temperature of the ferroelectric phase 
transition 7 . 



VIII. CONCLUSIONS 

The appearance of the spontaneous polarization in 
Sn2?2S6 compound is accompanied by the significant 
changes of electronic structure that are observed in all 
subbands of this crystal VB. At the transition from the 
paraelectric phase to the ferroelectric one, the significant 
changes also occur for the phonon spectra in the whole 
frequency range - for both external and internal vibra- 
tions of the crystal lattice. The complicate evolution of 
the energy spectra could be represented by following con- 
tributions. The fully symmetrical ('breathing') A g modes 
change the space distribution of the electron density in 
the elementary cell. This one initiates the stereochemical 
activity of the tin cations electron lone pair and support 
creation of their covalence chemical bonds with sulfur 
atoms (the hybridization of sp 2 d 5 type), and with phos- 
phorus atoms also. Thus in result of the ionic charges 
lowering, the short-range repulsion between tin and phos- 
phorus cations decreases, however the coulomb repulsion 
between tin cations still remains strong enough. These 
factors reflect the nonlinear interaction of A g B 2 type, 
they determine the anisotropy of polar deformations in 
the elementary cell and induce appearance of the dipole 
structure motives ( — Sn — (P2S6) — ), which are re- 
lated to the polar normal coordinates of B u symmetry. 
The dipole-dipole interaction correlates orientation of the 
local dipoles (pseudospins) and governs the spontaneous 
polarization of the crystal structure. In the structure re- 
arrangement, all low symmetry modes take part - the 
possibility of the nonlinear linking of A u B g B u type is 
obvious. This interaction correlates with involving of tin 
and phosphorus d orbitals into creation of the covalent 
bonds in Sn2?2S6 crystal lattice. Small difference of the 
paraelectric and ferroelectric phase's energies and acti- 
vation redistribution of the electron charge at the spon- 
taneous polarization appearance determine the presence 
of the three- well local potential in the Sn2?2S6 ferro- 
electrics. 
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